Layer Potential Techniques for the Narrow Escape Problem 

Habib Ammari* Kostis Kalimeris^ Hyeonbae Kang-^- Hyundae Lee^ 



Abstract 

The narrow escape problem consists of deriving the asymptotic expansion of the 
solution of a drift-diffusion equation with the Dirichlet boundary condition on a small 
absorbing part of the boundary and the Neumann boundary condition on the re- 
maining reflecting boundaries. Using layer potential techniques, we rigorously find 
high-order asymptotic expansions of such solutions. We explicitly show the nonlin- 
ear interaction of many small absorbing targets. Based on the asymptotic theory for 
eigenvalue problems developed in |2j , we also construct high-order asymptotic formulas 
for eigenvalues of the Laplace and the drifted Laplace operators for mixed boundary 
conditions on large and small pieces of the boundary. 
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1 Introduction 

An interesting problem in cellular and molecular biology is to estimate the mean sojourn 
time, also called the narrow escape time, of a Brownian particle in a bounded domain 17 
before it escapes through a small absorbing arc dQa on its boundary dQ. The remaining 
part of the boundary dQr = <9i7 \ dQa is assumed reflecting for the particle. The small arc 
often represents a small target on a cellular membrane. Its physiological role is to regulate 
flux, which carries a physiological signal [11], [1]. 

The narrow escape problem is connected to that of calculating the solution of a 
mixed Dirichlet-Neumann boundary value problem in $7, whose Dirichlet boundary is only 

a on the otherwise Neumann boundary. The escape time can be estimated asymptot- 
ically in the limit e := |(9f7a|/|(?^| — ?• by computing the asymptotic expansion of as 
e ^ 0. 

In this paper, we first consider a purely diffuse model. We provide mathematically 
rigorous derivations of the first- and second-order terms in the asymptotic expansion of 
the solution it^ as e — t- in the presence of a single or many small targets. When two 
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or more Dirichlet targets cluster together they interact nonhnearly. The clustering may 
affect significantly the asymptotics. Then we study the problem of eigenvalue changes due 
to the small targets. Finally, accounting for a drift term, we generalize our results to a 
mixed Dirichlet-Neumann boundary value problem for the drift-diffusion equation. 

The narrow escape problem of a free Brownian particle (without drift) through a 
small target was discussed in [5], [H], [15], and [16]. The method of [H], [15], [16] was 
generalized in [T2j to obtain the leading-order term of the solution to the corresponding 
mixed boundary value problem for the drift-diffusion equation. Matched asymptotics [17] . 
[Is] . [19] . [8] yield the expansion of the principle eigenvalue of the Laplace operator for 
mixed boundary conditions on large and small pieces of the boundary. The effect of 
clustering of the Dirichlet targets on the first eigenvalue was analyzed in [6] . The second- 
order term in its asymptotic expansion as the size of the target goes to zero was provided 
in [13] by determining the structure of the boundary singularity of the Neumann function 
for the Laplacian in a bounded smooth domain. In all of these papers, the derivations 
are quite formal. It is the purpose of this work to provide a rigorous framework for 
systematically deriving high-order asymptotic formulas for the solutions of the diffusion 
and drift-diffusion equations with mixed boundary conditions on large and small pieces 
of the boundary and the eigenvalues of the corresponding operators. Our derivations are 
based on layer potential techniques and the asymptotic theory for eigenvalue problems 
developed in [2]. 

While this work was in progress, we found out the existence of the work Some of 
results of this paper, especially those in Section 3.1, were also obtained in [9]. However, 
the method of derivation of asymptotics in [9] is formal and uses the method of matched 
asymptotics. 

This paper is organized as follows. In Section 2, we formulate the problem with and 
without the drift term, and review some facts on relevant Neumann functions and layer 
potentials on arcs. Section 3 is to derive (higher-order) asymptotic formula for the solutions 
to the narrow escape problem with a single or multiple (well-separated and closely located) 
absorbing regions. Section 4 is to deal with the eigenvalue perturbation problem in the 
presence of small absorbing region using the same layer potential techniques. Section 5 is 
for the narrow escape problem and the eigenvalue perturbation problem in the presence 
of the force field. The paper ends with a short discussion. 

2 Formulation of the problem 
2.1 Physical background 

Let be a bounded simply connected domain in M? with C^-smooth boundary. Suppose 
that dil. has two disjoint parts, the refiecting part dflr and the absorbing part dO,a, 
satisfying dO, = d^la U d^lr- Both dO,a and d^lr consist of finite number of open arcs. 
We assume that e := 190^1/2 is much smaller than 1 while \di}\ is of order 1. Here and 
throughout this paper \dQ\ denotes the arc- length of dQ. 

Suppose that a Brownian particle is confined to Q. The probability density function 
Pe{x, t) of finding the Brownian particle at location x at time t (prior to its escape) satisfies 
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the Fokker-Planck equation 



^{X,t) = Ape{x,t) - V ■ {F{x)pe{x,t)) 

with the initial condition 

Pe{x,0) = p{x), 
and the mixed boundary conditions for t > 

= on dQa, 

— peF • z/ = on dVtr- 

The force field F is given by F{x) = V0(x) for a smooth potential (p. The function p{x) 
is the initial probability density function; e.g., p{x) = for a uniform distribution or 

p{x) = 5y, a Dirac mass at y, when the particle is initially located at y. 

The function ^^(x) := Jq^°°|?£(x, t) dt, which is the mean time the particle spends at x 
before it escapes through dila, is the solution to 




' Av^-V ■ {Fv, 

Vs=0 
dVe 

du 



VeF-U = 



in Q, 
on dQa, 

on dQr- 



The function 'w^{x) 



V£{x)e '^^^^ is the solution of the adjoint problem 

-pe-'t' 



( Aws + F ■ Vw, 

We 

dw, 



e -r i ■ \ IMS 

'e 



du 







in il, 

on dO-a, 

on d^r- 



Suppose that p{x) = 6y. The function 



Ue{x) 



fnWe{x,y) dy 
e-'?^(f) dy 

is the solution to the mixed boundary value problem 



AUe + F • VUg 

Ue = 

dUe 

. dv 







in il, 
on d^ai 

on d^r- 



(2.1) 



The narrow escape problem consists of deriving the asymptotic expansion of as 
e — )• 0, from which one can estimate the sojourn time of the Brownian particle. A closely 
related problem is to construct the asymptotic of the principal eigenvalue of — A — F • V 
in with the mixed Dirichlet-Neumann boundary conditions for small e. 

The goal of this paper is threefold: (i) To find the asymptotic behavior of the solution 
Ufr to (|2.ip with F = (free Brownian motion) in the presence of one or multiple small 
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targets as e — )• 0; (ii) To find the asymptotics of the eigenvalues of —A in 17 with the 
Dirichlet boundary condition on dO,a and the Neumann boundary condition on dQr', (iii) 
To generahze these asymptotics to the case with nonzero drift {F ^ 0). 

We note that if the force field F = 0, then by integrating ()2.ip over Q the following 
compatibility condition holds: 



da 



(2.2) 



2.2 Neumann functions 

Let N(x, z) be the Neumann function for —A in corresponding to a Dirac mass at z G 0. 
That is, is the solution to 



/S.zN{x,z) = —Sx, x,z ^VL, 



dN 



1 



an 



N{x, z)da{z) =0, x G il. 



(2.3) 



The Neumann function is symmetric in its arguments; N{x,z) = N{z,x) for x ^ z £ 
It furthermore takes the form 



N{x,z) 



1 

'2^ 



ln\x — z\ + Rq{x, z), x,z£Q,, 



(2.4) 



where Rni', z) belongs to i/^/^(0), the standard Sobolev space of order 3/2, for any z G il, 
and solves 

A^Rnix, z) = 0, X G 17, 



dRn 



dfx 



(2.5) 



Suppose z G (917. Since dO. is smooth, only a half of any sufficiently small disk about 
z is contained in 17, and hence the singularity of N{x, z) is — (l/vr) In |x — Therefore, 
N{x, z) for z G (917, which we denote by Nq^, can be written as 



Nq^{x,z) = \n.\x — z\ + Rq^{x , z) , X G 17, z G (917, 

vr 



(2.6) 



where RQn{-,z) solves the problem 



AxRdn{x,z) =0, X G 17, 
dRdn 



dvx 



' xGai7. 



(2.7) 



|517| vr \x — z\ 



Note that the Neumann data in (|2.7p is bounded on 917 uniformly in z G 917 since 917 
is C^-smooth, and hence Rqq,{-,z) belongs to 77^/^(17) uniformly in z G 917. 
Let 

g{x):= I N{x,z)dz, x G 17. (2.8) 
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Then Ag = -1 in n and ^ = -\i}\/\dn\ on dn. Therefore one may use Green's formula 
and the mixed boundary conditions in (j2.ip with = to have 



u,{x) = g{x) + / iVan(x, z)^^da{z) + C,, x€ Q, (2.9) 
Jaua 

where 

n. = 

\dn\ 



In view of (12.61). (12.91) becomes 



Ce = j^i Ue{z)da{z). (2.10) 



ueix) = gix) - - [ ln\x - z\^^^da{z) + [ Rg^{x, z)^^^dc7iz) + C^. (2.11) 

If is the unit disk, then 

{x — z,Vx) {x — z,x) 1 — x-z 1 

\x — + |zp — 2x ■ z 2 — 2x ■ z 2 

for x E dQ. Therefore, ^gj'" \ xedn = 0, and hence Rgfi{x,z) = const. Because of the 
condition Jq^ N{x, z)da{z) = 0, we have Rgfi{x,z) = for all x & and z G dQ, and 
hence 

Ndn{x,z) = ln\x - z\, x £ Q, z e di}. (2.12) 

vr 

We also have 

g{x) = [ N{x, z)dz = -{I - |x|2), (2.13) 
which can be seen using Green's formula. In fact, we have 

N{x,z)A\z\'^ - A,N{x,z)\z\'^dz= [ (2N{x,z)\z\-^^^^^^\z\A da{z). 

Jan V (^^z J 

2.3 Single-layer potential on an arc 

Our aim is to derive an asymptotic expansion of du^{z)/di' on d^la for e small enough. 
We write this density function as the solution of an integral equation (|2.11|) on d^la- For 
doing so, we introduce two Hilbert spaces: For e > 0, define 

Xe:={^p:J Ve"^ - x^\^p{x)\'^ dx < +00} (2.14) 

with the norm 

ll^llx, = ( r Ve^-x^Hx)\'dxy/\ 



Then one immediately has 

||<^||xe = ll^llxi, (2.15) 
where tf(x) := eip{ex). We also define 

Y,:={i^GC'{[-e,e]):i^' GXe}, (2.16) 
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with the norm 

where ip' is the derivative of in the sense of distribution. 

We also recah the fohowing lemma; see [21 Chapter 5] and the references therein. 

Lemma 2.1. The integral operator C : Xi i— Yi defined by 

C[^]{x) = j \n\x-y\^{y)dy (2.17) 
is invertible. For a given function ip G Yi, G Xi is given by 



i:-\^\{x) = - / ^ lUiAllldy- ^7 ^ 2.18 

7r2 VI -x^ J-i x-y 7r(ln 2) VI - 



/or X G (—1, 1), where the constant a{ip) is defined by 

1 ^JY^^\z) 



a{i{j) = ip{x) + C 



vr 



2. 



Vi-y^ •''-1 2^ ~ ^ 



(x). (2.19) 



Note that, in view of ()2.18p and (I2.19p . we have 



C-\l]{x) = ——^==. (2.20) 



7r(ln2)\/l -x2 



3 Narrow escape of a free Brownian particle 

We now consider the narrow escape problem in the presence of a single or multiple (clus- 
tered or well-separated) absorbing arcs. 

3.1 A single small target 

Let x{t) : [—£, e] — s- be the arclength parametrization of dCla] namely, x is a C^-function 
satisfying |x'(t)| = 1 for all t G [— e,e] and 

dna = {x{t) ■.te[-e,e]}. (3.1) 

Since Ue = on d^la, it follows from ()2.1ip that 



/W-^/ ln\xit) - x{s)\^,{s)ds + j rit,s)ipeis)ds + Ce = 0, (3.2) 

where 

fit):=gix{t)), ^,(t) := r{t, s) := Ran{xit),xis)). 

By a simple change of variables, we obtain that 

— / ln\x{et) — x{es)\(p£{s)ds — / r{et,es)ifs{s)ds = f{et) + C£, (3.3) 
^J-i J-i 
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with ifeit) = eipe^et). 

Since the compatibility condition (|2.2p reads 



we may rewrite (13. 3p as 



1 

^e{s)ds = -\ni (3.4) 



-(£ + eCi)m{t) - + r(0, 0)|f]| = f{et) + C„ (3.5) 

TT TT 



where the operator Ci is given by 



. r n \xiet) - x{£s)\ , , , A , , , 

Moreover, since 

\x{et) - x{es)\ = e\t - s|(l + 0(e)), 

one can see that Ci : Xi i— )• Yi is bounded independently of e. Note that / is on Q 
and hence we may rewrite ()3.5p as 



(£ + e£i)[(^e](t) = |f7|lne-7rr(0,0)|Jl| + vrC^ + 7r/(0) + 0(e). (3.6) 
It then follows from ([2:20]) that 

1 



feit) 



'^'^''^ -r-(0,0)|f]|+C, + /(0) 



vr 



(In i) 



+ 0(8), (3.7) 



where the remainder 0{e) is with respect to the norm || • Hxi- 
Now, plugging ()3.7p into (I3.4p we arrive at 

(7, = M In - + r(0, 0)|1^| - /(O) + 0{e) (3.8) 
vr e 

and, consequently, 

^,(t) = -^-^=L= + 0(e), (3.9) 
where 0{e) is with respect to || • ||xi- Therefore, we obtain 

^,{t) = -\^—±= + 0{e), (3.10) 

where 0{e) is with respect to || • Hx.- 

Finally, combining (|2.1ip and (j3.10p yields 

Mx)=.(x)-H r ^:^^^^,,+c.+o(.) 

J-£ Ve^ - 
= g{x) - \n\N9u{x,x*) + Cs + 0{e) 

for X G provided that dist{x,dQa) > c for some constant c > 0. Thus we have the 
following theorem. 
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Theorem 3.1. Suppose that dQa is an arc of center x* and length 2e. Then the following 
asymptotic expansion of holds 

n,(x) = ^ln- + $n(2;,x*) + 0(e), (3.11) 
vr e 



where 



f N{x,z)dz-\Q.\NQn{x,x*)- [ Nix*,z)dz + \n\R9n{x*,x*). (3.12) 
Jn Jn 



The remainder 0(e) is uniform in x £ Q satisfying dist(x, (?ila) > c for some constant 
c > 0. Moreover, if x{t), —e<t<e, is the arclength parametrization of d^la, then 

, , ^^ \^\ 1 X 

where 0(e) is with respect to \\ ■ ||xj:- 

The leading order term in (jS.lip was obtained in [14] and the formula (13. IIP was 
obtained in [9j in a formal way (see (2.14) in that paper). We emphasize that the derivation 
of this paper is mathematically rigorous. Moreover, the method of derivation is quite 
different from that of [9]; It is based on the layer potential technique. We note that the 
corrector function ^q(x,x*) solves the following problem 

A^^n(x,x*) = -1 in 0, 
— = — |SZ|Oa;* on oil, 

which shows that (j3.1ip is nothing else than a dipole-type approximation. 

If $7 is the unit disk centered at 0, one can easily see from (|2.12p and (j2.13p that 

^q(x,x*) = ln|x - x*\ + ^(1 - |xp). 
Thus we have the following corollary. 

Corollary 3.2. Suppose that VL is the unit disk and dQa is an arc of center x* and length 
2e. Then the solution is given asymptotically as 

Ue(x) = In ^ + ^(1 - |xp) + In |x - x*| + 0(e), (3.14) 

where the remainder 0(e) is uniform inx £ Q satisfying dist(x, dQa) ^ c for some constant 



c. 



The formula (|3.14p was also obtained in |15] (and [9]). 
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3.2 Two well-separated small targets 

Let be a unit disk centered at 0. Suppose that dQa consists of two parts, say dfli and 
8^2, satisfying d^li n 80.2 = O- 

Let u be the solution to (12. ip with F = 0. In view of ()2.1ip . u is now written as 

u{x)=g{x)+ [ Nan{x,z)^^da{z)+ [ Nan{x, z)^^da{z) + C^. (3.15) 

Suppose that 8^li, i = 1,2, are parameterized by 

x{t) = (cost,sint), t e {si - Si, Si + £i). 
Then the (Dirichlet) boundary conditions on 8rta yield the following integral equations 

0= / N9n{x{s),xit))cP{t)dt+ / N9n{x{s),x{t))(j){t)dt + Ce (3.16) 

-'Sl— £l J S2—e2 

for s £ {si- Ei, Si + Ei) where (j){t) := ^{x{t)). 

Recall that Nqq{x,z) = — ^ In \x — z\, x G 80, z G 8Q for the unit disk. If we make 
changes of variables s i— )■ Sj + e^s and t ^ Si + Eit for s,t £ (sj — Ei, Si + Ei), the integral 
equations (|3.16p become 

f 

/ In \x{si + Eis) — x{si + Eit)\ ipi{t) dt 

+ y In |x(si + Eis) - x{s2 + £201 'f2{t) dt = ttCe, 
J In \x{s2 + E2s) - x{si + Eit)\ipi{t) dt 

+ J In \x{s2 + E2S) - x{s2 + £2^)1 ¥^2(0 dt = ttCs, 

where (/7j(t) = Ei4>{si + Ejt) for t G (—1, 1). 
Let 

1 

Ci=- ipi(t)dt, i = l,2. (3.18) 

Because of (|2.2p . we have 

Ci + C2 = -1. (3.19) 

One can easily check that 

In \x{si + Eis) - x{si + Eit)\ = lnEi\s - t\ + 0{Ei), 

( \ (3-20) 

\n\x{si^E^s)-x{^Sj^Ejt)\ = In \x{si) - x{s2)\ + O ( j , i^j. 

Therefore, the system (|3.17p can be written in the form 

{A + B)(^^=.(''^- - ^^1^(7) - f , (3.21) 

' \^2j \Ce-Ci\n\x{si) -x{s2)\-C2\nE2j 



(3.17) 
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where the operator A : Xi x Xi i— t- Yi x Yi is given by 



and 5:XxXi->YxY satisfies 

/ 

B = 

' O 



A:-- 



Oiei) 



C ' 



o 



y \^|x(si)-a;( 



S2j 



\x{si)-x{s2)\ 

0{e2) 



in the operator norm. Here C is given by ()2.17p . 

Let us now assume that \x{si) — x{s2)\ > c for some c > 0, in other words, the two 
targets are well-separated. It then follows from (j2.20p that 



Ce - Cilnei - C2ln|a;(si) - a;(s2)| / 1 , r^r i 2 , -i 



Ini 



\VT~- 

- Ci In |x(si) - x(s2)| - C2 lne2 / 1 



In- 



+ 0{Jei + e 



(3.22) 



2J ' 



where ©(yef + ef) is with respect to || • ||xi- 

Integrating (|3.22p over (—1,1), we obtain a system of two equations for Ci, C2 and 



Ci In ^ + C2 In \x{si) - xis2)\ - = 0(^/^+4), 
Ci In |x(si) - x{s2)\ + C2 In ^ - C7, = 0{^/4 + 4), 



which together with (j3.19p yields 



Thus we have 



£2 



Ci 
C2 

Ce 



2|x(si)-a:(s2)| 
In £l£2 

4|x(si)-x{s2)P 



+ OUei + ei), 



El 



2|x(si)-a-(s2)| 
In £l£2 

4|x(si)-x(s2)P 



+ o(V^f + e^), 



(ln|x(si) -x(s2)|)2 -Inf Inf 

In ^1^2 

4|x'(si)-a'(s2)P 



9'i(s) 
9^2(5) 



In 



£2 



2|x(si)-x(s2)| 



1 



In 
In 



£l£2 



4|a;(si)-a;(s2)|^ Vi - S 
£1 

2|x(si)-x(s2)| 



1 



+ 0{Jei + ei), 



1„ £1£2 

4|x(si)-a;(s2)P 



+ 0(Ve^ + e|), 



(3.23) 
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where ©(yef + ef) is with respect to || • ||xi- By scahng back, the following expansion of 
(j){s) holds as e — )• 0: 



In 



£2 



In 



2\x{si)-x(s2)\ - , /^r / 2 , 2\ ^ I I ^ 



£l£2 l~2 1 



In .,^.„ f^.„ . 1 



By substituting this formula to p.l5p . we obtain the following theorem. 

Theorem 3.3. Let dO,i, i = 1,2 be parameterized by 

x{t) = (cos t, sinf), t £ {si — Ei, Si + e,). 

Assume that there is a constant cq > such that \x{si) — x{s2)\ > cq. Then the solution 
u to h2.1\) is given asymptotically by 



(ln|x(.i)-x(.2)|)^-lnf Inf ^ 1 
= 1^ £i£2 + - 1^1 ) 

In . 1 , . .1 In ^ 



+ . In \x - x{s,)\ + In |x - x(s2)| + el + e^), 

i\x{,Sl)-x{s2)\^ i\x(si)-x{S2)\' 

(3.24) 



where 0(\/ + £2) ^'^ uniform in x £ Q satisfying dist(x, 90iU9r22) ^ c /or some constant 
00. 

If e := ei = £2) then (j3.24p takes the following simpler form: 

1,2 1, , ,9, 1, |x — |x — x(s2)| ^, , , , 

u{x) = - In - + - 1 - X 2 + In ' I ; V r ^ + ^ ^ • 3-2^ 

Formula (|3.25p shows that the leading order term of the escape time when there are two 
targets of the same size is reduced to half of that when there is a single target, which is 
quite natural. The second-order term (dipole type) shows that the escape time is reduced 
to half of the sum of two dipoles increased by the relative position of the targets. 

3.3 Two clustered small targets on the unit circle 

We now consider the case when there are two clustered targets. We suppose that 

• ei = e2(=: e), 

• S2 — si = de, d > 2. 

We still assume that is a unit disk centered at 0. 
Let 



Ci[(l)]{s) := J In \d + s - t\(p{t)dt, 
C2[(t>]{s) := j \n\d- s + t\(t){t)dt. 
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The system of equations (I3.17P now takes the form 



where A, i3:XxXi— t-YxY are given by 

A: 



C £2 
£1 C 

and B = 0{e) in the operator norm. 

The foUowing invertibihty result follows from [101 



Lemma 3.4. The operator ^:XxXi— s-YxY is invertible. 
It follows from Lemma 13.41 and ()3.19p that 



= AC, + Ine) I^A-' (I) + 0(e)) . (3.27) 
Define aj = aj{d), j = 1, 2, by 

Then, we have ai = a2(= a)- In fact, if {ipi,(p2)'^, where T denotes the transpose, is the 
solution of 

(::)-(;)■ 

so is {ip2{—t),ipi{—t))'^, and hence we have 

Mt) = ^2{-t). (3.30) 
Integrating (|3.27p over d^la yields 

which combined with (I3.15P and (|3.27|) . gives the following result. 
Theorem 3.5. Let dfli, i = l,2 be parameterized by 

x(t) = (cost, sin t), t £ (si — £, Si + e) 
with \si — S2I = de. Then the solution u to ^2.1\) is given asymptotically by 



u{x) = In ^ - + ^(1 - \x\^) + ^ In |x - x{si)\ + ^ In |x - x(s2)| + 0(e), (3.31) 

where 0(e) is uniform in x G il. satisfying dist(a;, U 91^2) ^ c for some c > 0. 
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In view of (f330]l . it follows from (f3:29]l that 



or 



j [ln\s-t\-ln\d-s-t\](pi{t)dt = l. 



By taking s = 1 one can see that a{d) — )• oo as d — )• 2+, and hence the solution converges 
to the one with an absorbing arc of length 4e (modulo 0{e)). 

It is interesting to compare the formula ()3.3ip with ()3.25p . The nonlinear interaction 
between the two small targets is described by the term In ^ — 2a(d) ™ (!3-3ip while it is 
described by the term —^\n\x{si) — x{s2)\ in (j3.25p . 

3.4 Multiple small targets on the unit circle 

We now extend the above analysis to larger clusters. We consider clusters dili, i = 1, . . . ,n, 
parameterized by 

x{t) = (cos t, sint), t £ {si — e, Si + e). 

We assume that each dili is separated from its neighbors by a distance comparable to e. 
Let dij := .^i— ^. Define the operators Cij and A by 



and 

Let 



Cij[4>]{s) := j \n\dij + s - t\(j){t)dt 

A . — {^ij^nxn- 













\-tr 




dt 











Following the same lines as in the proof of Theorem 13.5^ we can easily prove the 
following result. 

Theorem 3.6. Let dO,i, i = 1, . . . ,n be parameterized by 

x(t) = (cost, sin t), t £ {si — £, Si + e) 
with Si — Sj = dij£. Then the solution to \2. 1\) is given asymptotically by 

1 1 " 

u{x) = - Ine + - x{si)) + 0(e), (3.32) 

ai H h an n -f^ 

where the remainder 0(e) is uniform in x £ with dist(x, U"^^9rij) > c for some constant 
c> 0. 
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Note that 

1 If 1 

-y^,^n{x,x{si)) ^ vq{x) := —— ^n{x,z)da{z) as n = O(-) +00, 
n ^ \on\ Jq^ £ 

and the corrector u q satisfies the equation 

AfQ = — 1 in ri, 

with, as already known in the homogenization theory [3], the effective Neumann boundary 
condition 

4 Asymptotics of the eigenvalues 

Let = Ag^^ < Aq'^^ < Aq'^^ < ... be the eigenvalues of —A on with Neumann conditions 
on dQ. Let Uq"'^ denote the normalized eigenfunction associated with Ag"'^; that is, it 



satisfies l^^c/^P = 1- Let a; ^ {y AQ''^}j>i. Introduce N^{x,z) as the Neumann function 
for A + w'^ in $7 corresponding to a Dirac mass at z. That is, is the unique solution to 



(A^. +a;2)iV'"(x,z) = -(5^ in 0, 

dN^ (4.1) 
— — = on oil. 

ov 

We recall two useful facts on the Neumann function [2, (2.28) & Lemma 5.1]. 
Lemma 4.1. The following spectral decomposition holds pointwise: 

4^\x)ui^\z) 



iV-(x,z) = j;^^4P^, x^zen. (4.2) 
j=i ^0 ~^ 

Moreover, if Xq°^ is simple and V is a complex neighborhood of \J X^^^ which contains no 

•^/a^ other than \J Xq°\ then for 00 G V the Neumann function on the boundary d^l, 
which we denote by Nq^, has the form 

N^n{x,z) = --\n\x-z\+ '^J)'J ^ +R%n{x,z) forxen,zedn, (4.3) 

where RgQ{-,z) belongs to H^^'^{Q) for any z € dQ and uj 1— )■ Rqq^ is holomorphic in a 
complex neighborhood of Xq°^ . 

Suppose that Aq'"'^ is simple. Then there exists a simple eigenvalue A^"'"^ near Aq''"^ 

associated to the normalized eigenfunction ui''°^ {J^ = 1) satisfying the eigenvalue 

problem 



-Au^i"^ = Xi^'^i'"^ in n, 
4^0) = on dQa, 

— - — = on dQr- 

< ov 
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(4.4) 



Our aim in this section is to derive a high-order asymptotic expansion 
as e — ^ 0. 

Set = dui''°^ /du on dQ,a- By Green's formula, one can easily see that 

Ueix) = {X'^J°'^ - uj"^) [ 4^°\z)N'^{x,z)dz+ [ ipeiz)N'^ix,z)da{z), 

Jq JdUa 

which shows that 



iPe{z)Ny^^"'\x,z)d(jiz) = on d^a- 

IdQa 

Therefore, the eigenvalue problem (j4.4p reduces to the study of the characteristic values 
of the operator- valued function A i— )• ^^(A) given by 

AeiXMixit)) = I' ^{s)N^{x{t),xis))ds, 

where x{t) is defined by (j3.ip . 

Let y be a neighborhood of Aq''"^ in the complex plane such that Aq''"^ is the only 
eigenvalue of —A in $7 with Neumann boundary condition. The following pole-pencil 
decomposition of ^e(A) : — Yg holds for any X £ V: 

AeiX) = --Cs + -^^+n,iX), 

^ X}^^"' - X 
where ^ 

Ce[^p]{x{t)) = J In \s — t\ip{s) ds, 

fCe is the one-dimensional operator given by 

K,M{<t)) = ui%x{s)Mx{s))ds^ <{xit)), 

and 

7^e(A)M(x(^)) = £ Rg{x{t),x{s))^{s)ds. 

An application of the generalized Rouche theorem [2] shows that the operator-valued 
function Ae{X) has exactly one characteristic value in V. Theorem 5.7 in [2] gives its 
asymptotic expansion for e small enough. The following theorem holds. 

Theorem 4.2. We have 

xOo)_xi^o) ^ f C-\uf\x{t))u^i^^\x{t))dt 

h (4-5) 

+vr2y ^£-l7^,(A(^"))C^4'"^](^(^))4'"^(^(^))^^■ 

In particular, the leading-order term in the expansion o/ A^"'"'' — Ag''"'' is given by 

aOo)_aOo)«__I_|4.o)(^*)|2^ (4.6) 

where x* is the center of dQa- 

Note that if Aq'"^ is a multiple eigenvalue of —A in $7 with Neumann boundary con- 
dition, then in exactly the same way as Theorem 3.24 in [2j, we can address the splitting 
problem in the evaluation of eigenvalues of (j4.4p . 
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5 Narrow escape problem in a field of force 

We now deal with the case when the force field F = is not zero. As before, we 
consider the problem of estimating the escape time and that of deriving asymptotics of 
the perturbed eigenvalue. 

5.1 Asymptotics of the solution 

Let the Neumann function N^{x, z) be defined by 



A,N^{x, z)-V,- {F{z)N^{x, z)) = -6^, x,zen, 



F ■ u,N^ 



jean 



\dn\' 



(5.1) 



(5.2) 



Ue{x)=[ N^{x,z)dz+ [ ips{z)Ng^{x,z)da{z) + [ Ue{z)da{z), x e 

where ipe = du^jdv on d^a- Therefore, the density function satisfies the integral 
equation 

/ In |x — z|99e(2;) (i(T(z) + / R^Q{x^z)ip^{z) da{z^ = — / N^{x,z)dz — C£ 

for X £ dQa, where 



The (boundary) Neumann function, denoted by Nq^{x,z), has the form 

Nqq{x,z) = In \x — z\ + RgQ{x, z), X G ri, z G (9il, 

vr 

where RgQ{-,z) belongs to H^^'^{Q) for any z G di}. 

The solution to (j2.ip admits the following integral representation: 



a 



1 



Us{z) da{z). 



Moreover, since F = Vcp, we have 

V • (e'^Vue) = -e'l' in fi, 
and hence ip^ satisfies the compatibility condition 

e'^("Ve(2)d^(^) = - / e-^^^^dz. 
an Jn 

In exactly the same way as Theorem 13. 11 we can prove that 

Ine 

In V ^ 



^)dz[ — -RUx\x*: 



+Ce+ [ N^ix*,z)dz 

Jn 



1 



(lnl)^/r^ 



+ 0(8), 



where x{t) is defined by ()3.ip and x* is the center of dila- Consequently, we have 
Ce=(^j e'^^''^-*^''*Uz^(^^ln^ + R^^{x*,x*)^ - j {x* , z) dz + 0{e), 
which yields the following result. 
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Theorem 5.1. The following asymptotic formula 

u,{x) = l(^J e'^^^)-'^^^*)^ dz\n^ + X*) + 0{e), 

holds uniformly for x E CI satisfying dist(a;, dQa) > c for some constant c > 0. Here 
is given by 



<P^{x,x*) = [ N^{x,z)dz- [ N^ix* 
Jn Jn 



X* , z)dz 

+ ( j^e't'^''^-t'^^*Uz]{R^9^{x\x*)-NUx,x*)). 



(5.3) 



5.2 Asymptotics of the eigenvalues 

Consider the eigenvalue problem 



Ane,F + F ■ VUe,F = -^e,FUe,F, 
Ue,F = on dQ,ai 

^ =0 on oilr- 
ov 



in Jl 



(5.4) 



Since F = V<j) for a smooth potential function we get 

e-'^V • {e^Vu) = Au + F-Vu, 

and hence we can realize the Neumann Laplacian with a drift -F in as a non-negative, 
self-adjoint operator on L^(0, e^'^dx) via the closure of the Priedrichs extension of the 
nonnegative quadratic form 



(u,v) / Vu-Vve'^dx. 
Jci 



Let = X^p < Aq^];, < Aq^];, < ... be the eigenvalues of — A — F • V on Q with Neumann 
conditions on 50 and let Wq^];,, u^p, u^p, ... be the corresponding eigenfunctions with 



Id I'^o'fI^^'^ ~ ^- Define for a; ^ '^o'''}j>i the Neumann function by 



A.ATf (x, z) - V, • (F(zX (x, z)) + iv^N^ (x, z) = -5,,, x, z G Q 



F- u,N^ 



0. 



(5.5) 



Set 



From 



</^£,F := — ^ on dQ.a. 



/ ipe,F{z)N^r- — {x,z)da{z)=Q and^a, 
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it follows that the eigenvalue problem ()5.4p reduces to the study of the characteristic values 
of the operator- valued function A i— )• Ae^ri^) given by 

AsA>^Mi^it)) = f v{s)N^{x{t),x{s))ds. 

Following the same arguments as for the case without drift, we can prove the following 
result. 

Theorem 5.2. // Aq"*^^ is simple then there exists a simple eigenvalue X^Jp of con- 
verging to Xq"p as e goes to zero and the following formula holds: 

ti->^tF^-^yo:X-*)?e^^^*'- (5.6) 

6 Conclusion 

In this paper, we have provided mathematically rigorous derivations of the first- and 
second-order terms in the asymptotic expansion of the solutions of diffusion equations in 
the presence of a single or many small targets. We have shown the nonlinear interaction 
between many targets. We have also studied the problem of eigenvalue changes due to 
the small targets using the approach developed in [2j. The generalization of the results 
of this work to nonsmooth domains containing corners and cusps will be the subject of a 
forthcoming paper. 
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